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B s -B s Mixing in Supersymmetric Grand Unified Models 
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We study B 3 -B s mixing in Supersymmetry grand unified SO(10), SU(5) models where the mixings 
among the second and third generation squarks arise due to the existence of flavor violating sources 
in the Dirac and Majorana couplings which are responsible for neutrino mixings. We find that when 
the branching ratio of r — > /17 decay is enhanced to be around the current experimental bound, 
B a -B a mixing may also contain large contribution from supersymmetry in the SO(10) boundary 
condition. Consequently, the phase of B s -B s mixing is large (especially for small tan f3 and large 
scalar mass mo) and can be tested by measuring CP asymmetries of B s decay modes. 
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Flavor changing processes are important not only to 
test the Kobayashi-Maskawa theory |l[ and to determine 
the parameters in CKM (Cabibbo-Kobayashi-Maskawa) 
matrix, but also to examine new physics. Recent mea- 
surement of B s -B s mass difference, 



AM S = 17.77 ± 0.10 ± 0.07 ps" 



(1) 



by D0 and CDF collaborations [2| can impact on the fla- 
vor structure of new physics beyond the standard model 
(SM) 0,i. The experimental constraints for new physics 
are not very severe yet since deviations from the SM pre- 
diction can be buried in the errors of CKM parameters 
and lattice calculation. In other words, there is still room 
for new physics. However, the parameters are expected 
to be determined more accurately in the near future. Be- 
sides, the CP asymmetry of B s decay can be observed by 
direct measurements, and we will get important informa- 
tion of new physics. 

Supersymmetry (SUSY) is one of the most promising 
candidates of new physics. SUSY can provide a natural 
prospect to have a large hierarchy in the theories, and 
in the minimal SUSY standard model (MSSM), gauge 
forces can unify at a high scale, which leads to a success- 
ful realization of grand unified theories (GUTs). How- 
ever, the flavor sector has not yet been well accepted in 
the MSSM due to the fact that SUSY breaking terms 
can induce large flavor changing neutral currents (FC- 
NCs) . Actually, the experimental constraints of FCNCs 
introduce flavor degeneracy of the SUSY particles, espe- 
cially in first and second generations, if SUSY particles 
are lighter than around 2-3 TeV 

In order to suppress the SUSY FCNCs, the squarks 
and sleptons are assumed to be degenerate at the GUT 
scale as a nature of the SUSY breaking. Even though 
the degeneracy is realized at a scale, the renormalization 
group equation (RGE) flow induces flavor violation for 
squarks and sleptons at low energy. In this scenario, 
the flavor violating effects in SUSY breaking terms are 
small at the weak scale and satisfy the current FCNC 
constraints. The small flavor violations originate from 
mixings in the Yukawa couplings characterized by the 



CKM mixings as well as the neutrino mixings. In the 
MSSM, the induced FCNCs in the quark sector are not 
large since the CKM mixings are small. In the lepton 
sector, on the other hand, sizable FCNC effects can be 
generated and a testable amount of flavor violating lepton 
decay can be obtained Q, which is related to the large 
mixings for the neutrino oscillations. 

In the grand unified models, the flavor violation at the 
weak scale can be related to the GUT scale physics. Ac- 
tually, as a consequence of the quark-lepton unification, 
the large neutrino mixings not only introduces flavor vi- 
olations in the lepton sector, but also in the quark sector 
as well. The relation of flavor violation in the quark and 
the lepton sectors depends on unification of matters and 
how to obtain light neutrino masses. Therefore, investi- 
gating the FCNC effects, we may obtain a footprint of 
the GUT models. 

The new result on B a -B s mixing can restrict the flavor 
violation in the quark sector involving b and s quarks. In 
grand unified models, the B s -B s oscillation can be cor- 
related to the t — > fjrf decay. Since Br(r — ► fij) is being 
measured at the B factories, the future results will be 
able to probe new contributions from the GUT models. 
In this letter, we calculate B s -B s mixing and Br(r — > (i<y) 
in SU(5) and SO(10) GUT models, and_study the im- 
plication of the correlation between B s -B s mixing and 
Br(r — > fij) decay in the CP asymmetry of B s -B s mixing 
(which is under experimental investigation) to decipher 
GUT models. 

The existence of second-third generation (23) mixing 
elements in squarks and slepton mass matrices generate 
B s mixing and r — > /17 decay respectively. We first in- 
vestigate the arise of 23 elements in squark and slepton 
mass matrices in the grand unified theories. 

SU(5) model: In a SU(5) grand unified model, the 
superpotential involving the Yukawa couplings is as fol- 
lows: Wy = Y^/AWaOjHs +_V2Y dlJ W t 5jH- 5 + 
YvijSiNjH*, + M vi j jl NiNj, where 5 contains the right- 
handed down-type quarks (D c ) and left-handed lepton 
doublets (L), and i, j denote the generation indices. The 
left-handed quark doublets (Q), right-handed up- type 
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quarks (U c ), and right-handed charged- leptons (E c ) are 
unified in 10 multiplet, and N is the right-handed neu- 
trino. Since two large neutrino mixings have been ob- 
served in nature, the Y v coupling is expected to have large 
off-diagonal elements, which will generate off-diagonal 
terms in the SUSY breaking scalar mass matrices for 
the scalar 5 multiplet via H§ and N loops @. One, 
thus, expects a large 23 element in the SUSY breaking 
scalar mass matrices for D c and L due to the large at- 
mospheric mixing and possibly large v T Dirac Yukawa 
coupling. Therefore, the SUSY contribution for the am- 
plitude of B s -B s mixing can be enhanced along with the 
branching ratio of r — > fi"f. 

Since our purpose is to investigate the flavor violation 
in the 23 sector of squark and sleptons and to probe how 
it relates the B s -B s mixing and flavor violating r de- 
cay, we are only discussing the 23 element and its effects 
in quark and lepton flavor violation in this letter. For 
this purpose, we consider the following simplified SU(5) 
boundary condition for the SUSY breaking scalar mass 
matrices at the GUT scale where SU(5) is broken to the 
SM: 
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The 5 denotes the flavor mixing term arising from the 
neutrino Yukawa couplings discussed above. We as- 
sume the above boundary condition in the basis where 
the down-type quark Yukawa matrix is diagonal (at the 
GUT scale): Y d = Y d diag , Y a = Vg KM Y* ia *P u V uR ,Y e = 
V e LY^ lag P e V^ R , where the up- and down-type quarks 
and the charged-lepton Yukawa couplings Y u ^ are real 
(positive) diagonal matrices and P u ,e are diagonal phase 
matrices. In a minimal SU(5) GUT, in which only H$ 
and H$ couple to matter fields, we have V u r — Vckm, 
V eL = V eR = 1, and Y diag = Y e dia s. We do not assume 
the minimal choice of Higgs fields, but assume V e L — 1 
to keep the relation of flavor violation between the quark 
and the lepton sectors. 

We note that one can also have first-second and first- 
third generations mixings in squarks and sleptons. How- 
ever, they can be small by the choice of 13 mixing in Y v 
and the hierarchical pattern of the neutrino coupling. If 
these elements are large, fi — ► ej will become large. 

SO(10) model: In a SO(10) model, the flavor vio- 
lations can be more enhanced compared to the SU(5) 
case since all matters are unified in spinor representa- 
tion and couple to 10 and 126 Higgs fields Q, e.g., 
W Y = |/iijl6il6j-10+ |/yl6il6jl26. The mixings in 
the neutrino Dirac Yukawa coupling may be small since 
right-handed neutrinos are also unified with other quarks 
and leptons. In a SO(10) model, however, there could be 
sources for large flavor violations in Majorana couplings 



for both left- and right-handed neutrinos in the type II 
seesaw scenario Q. The Majorana couplings are unified 
to the 16 • 16 • 126 coupling, and also affect the quark 
fields. The couplings will give rise to observable amount 
of flavor violations to the sparticle mass matrices via the 
GUT particle loops. Based on the above discussions, the 
following mass terms can arise in a SO(10) model at the 
GUT scale: 
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1 
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When the matters couple to only 10 and 126 
Higgs fields, the Yukawa matrices are symmetric and 
the boundary condition is Yd = Y dlas _P,i, Y„, — 



V L Y a ' 



l PuV CK M,Y e = V ql Y L dia s PeV T_ We do nQt ag _ 
sume the minimal choice of Higgs fields but assume 
Vqi ~ 1. It needs to be noted that the diagonal phase ma- 
trices P u ,d,e can not be rotated away for this boundary 
condition and these parameters enter into our calcula- 
tions which we will see later. 

It appears that in order to suppress the proton de- 
cay and to obtain the correct fit to fermion masses one 
needs to extend the above minimal SO(10) model. The 
new superpotential includes 120 Higgs field: Wy = 
^hijWiWjlO + 1/^16^16^126 + i/i^ 16,16,, 120. In 
this case, the symmetric nature of the Yukawa matri- 
ces is lost since 120 Higgs coupling h' is antisymmetric. 
In this context, the Hermitian Yukawa matrices can be 
considered with 120 Higgs and a parity symmetry to re- 
duce the number of parameters and to solve SUSY CP 
problem 0. The SO(10) symmetry is broken in the ba- 
sis where Yd is diagonal, since the left- and right-handed 
fields are rotated by conjugated unitary matrices, e.g. 
Q^VQ and U c -> V*U C . In the original SO(10) basis, 
the SUSY breaking mass matrices are real by the parity 
symmetry. As a result, in the basis where the down- type 
quark Yukawa matrix is diagonal, the squarks and slep- 
ton masses at the GUT scale are related by the following 
relation 
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(5) 



So we see that two mass relations are possible in a SO(10) 
model at the GUT scale. We call one of them symmet- 
ric (since the Yukawa couplings are symmetric) and the 
other one hermitian (since the Yukawa couplings are her- 
mitian). 

We use the above mass matrices for the boundary con- 
dition at the GUT scale and then calculate the masses 
at the weak scale by using RGEs to calculate the mix- 
ing of B s -B s and Br(r — > fj/y). In calculating the mass 
differences of mesons, one encounters non-perturbative 
factors originating from strong interaction. In the ratio 
of B s and Bd mass differences, many common factors 
cancel, and the ratio can be calculated more accurately 



3 



rather than the respective mass differences. We have 



AM S SM 
AMf M 



id 



(6) 



where £ = ^B^Jb s / {^B^J Bd ) = 1-23 ± 0.06 is a 
ratio of decay constants /s s(d) and bag parameters Bb sW 
for i? s (d) mesons. 

It is convenient to parameterize the SUSY contribu- 
tion by two real parameters C Rs and <ps, m model- 
independent way as 



C B ,e 



2,<i 



M 12 (B S ) 
M 12 (S S )SM 



(7) 



where M 12 (B S ) = M 12 {B S ) SM + M 12 (B s f VSY denotes 
the off-diagonal element of B s -B s mass matrix. Super- 
script SM (SUSY) stands for SM (SUSY) contribution. 
The mass difference is given as AM S = 2|Mi 2 (-B s )|- 

We now discuss the SUSY contributions in B s -B s mix- 
ing. When the flavor degeneracy is assumed in the 
MSSM, the chargino diagram dominates the SUSY con- 
tributions of Mi2(B 8i d). In this case, <j>s s — in Eq.©, 
and the ratio of mass differences in the MSSM is almost 
same as in the SM. In the general parameter space for the 
soft SUSY breaking terms, the gluino box diagram domi- 
nates the SUSY contribution of Mi 2 (B s ). The gluino (g) 
contribution can be written naively in the following mass 
insertion form 



Mf 2 



a l(5 d LL )l 2 + (6 d RR ) 2 32 ] - b (S d LL ) 32 (5 RR ) 32 , (8) 



where a and b depend on squark and gluino masses, 
and Sf LRR = (Mj)LL,RR/rh 2 (m is an averaged squark 
mass). The matrix Mj is a down- type squark mass ma- 
trix (Q, £) c t )Mj(Qt, D C ) T in the basis where down-type 
quark mass matrix is real (positive) diagonal. When 
squark and gluino masses are less than 1 TeV, a ~ 0(1) 
and b ~ 0(100). We also have contributions from 
8^ R , but we neglect them since they are suppressed by 
(mb/msusy) 2 - It is worth noting that the SO(10) bound- 
ary condition gives much larger SUSY contribution of 
Mi 2 (B s ) compared to the SU(5) case since both off- 
diagonal elements for LL and RR are large and b 3> a in 
the formula, Eq.®. 

For the calculation of observables at weak scale, we 
need to use the basis where Yd >e are real (positive) diago- 
nal matrix. In this basis, the boundary condition is given 
as Afg = P d M\P\, M\ ~ P e M\P\. Since the decay 
width is proportional to the squared absolute values of 
decay amplitudes, the phase of 6 in the SUSY breaking 
mass matrix at the GUT scale and the phase of P e in the 
Yukawa couplings are less important for Br(r — > /17). 
On the other hand, the phases of S and Pd are important 
for Mi 2 (B s ) slJSY . Due to those phases, the argument of 
Mi 2 (B s ) slJSY can be completely free. 
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FIG. 1: Minimal and maximal values for ratio of mass differ- 
ences versus Br(r — > y/y) under SU(5) and SO(10) boundary 
conditions. We show |5| = 0.05,0.1,0.15 points. The dotted 
lines show 90% CL region of the experimental data. 



In Fig.l, we plot the maximal and minimal values of 
the ratio of mass differences versus branching ratio of 
t — ► fj/y in the case of tan/3 = 10 where tan/3 is a ratio 
of up- and down- type Higgs vacuum expectation values. 
In the plot, we use Mxi 2 = 300 GeV and A = for the 
universal gaugino mass and the universal trilinear scalar 
coupling coefficient at GUT scale. The universal scalar 
mass at the GUT scale is mo = 200 for the SU(5) plot, 
and m = 200,400 GeV for the SO(10) plots. We use 
\Vtd/V t s\ = 0.192±0.009 which is obtained by the global 
CKM parameter fit without using experimental data for 
AM S [llj . The ratio of mass differences is proportional 
to £ 2 /|V t d| 2 . The Br(r — > /17) is almost proportional 
to tan 2 f3, while AM S / AMd does not depend on tan/3 
much. We find that Br(r — > ^7) does not have much 
dependence on tuq (for mo < 500 GeV), while the SUSY 
contribution of AM S depends on mo- 

In order to illustrate that any phase is possible for 
Mi 2 (B s ) slJSY , we plot the real and the imaginary part 
of M\ 2 {B S ) in the case where Br(r — * fjrf) = 6.8 x 10~ 8 
[12} for tan/3 = 10 in Fig. 2. We use the same values for 
M1/2 and A as in Fig.l and m = 500 GeV. Using SU(5) 
and SO(10) (hermitian) boundary mass values, we vary 
the phase of S. With SO(10) symmetric boundary condi- 
tions, we fix S to be real (positive) and vary the phase in 
the diagonal phase matrix P d - In the case of SO(10) with 
hermitian boundary conditions, the plot is a double-circle 
due to the gluino contribution as shown in Eq. {8} . Due to 
the RGE effect, the double-circle does not overlap com- 
pletely. We note that even if the radius of circle becomes 
large, AM S — 2|Mi2(B s )| has experimentally allowed 
solutions, as long as |Mi 2 (B s ) SUSY | < 2|M 12 (B S ) SM |, 
though one needs to adjust the phases in boundary condi- 
tions. It is worth emphasizing that the phases of Mi 2 (B s ) 
are large in such solutions. 
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FIG. 2: Re-Im plot for 2Mi2(B s ) when Br(r — > firf) saturates 
experimental bound. We use yj Bb s /s s = 262 MeV [13]. The 
dotted lines show 90% CL region of the experimental data. 
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FIG. 3: max(2 |</>bJ) (degree) versus Br(r — » fj/y) 



In order to see the maximum allowed phase in the 
B s -B s mixing, we plot the maximal value of 2|</>bJ 
versus Br(r — > ^7) in Fig. 3. We use the same val- 



ues for M1/2 and A 



Fig.l. 
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the boundary condition becomes large, the radius of 
the circle in Fig. 2 becomes large, and \4> Bs \ can be 
large consequently. One can approximately obtain 
max(sin2 \<j> B .\) - \M? 2 /M^\. When the \<j> Bs \ is max- 
imized, we find C Bg — cosmax(2 \4>b 3 |)- The model 
independent constraints for C Bs and <j) Ba are C Bg — 
0.97±0.27and 2cj> Bs = (-4±30)°U (186±30)° p). The 
phase 4> Bs can be measured by CP asymmetry of the de- 
cay B s — ► J/ijj(t> and the semi-leptonic decay B s — > 
The phase ^>s a in the SO (10) case can be larger than 
in the SU(5) case since both LL and RR elements can 
be large in the SO(10) case. As depicted in Fig. 3, there 
is a chance for 2<f) Bg to be around 20° in the SO(10) 
case before the parameter space gets ruled out by the 



Br(r — > /17). In order to obtain a large phase, large mo 
and small tan f3 are needed, which leads to an important 
implication. The Higgs mass bounds for MSSM restricts 
the lower values of tan/? and Mi/%. In the minimal su- 
pergravity model, the scalar mass mo is restricted to be 
less than around 200 GeV (for m < 1 TeV) [H for 
tan/3 = 10 by the WMAP data Q, and as a result, the 
phase \4>b s I can not be very large. Interestingly, the large 
mo solution (mo > 1 TeV) for dark matter content [l5j 
may generate large 4> Bs ~ 90° which is allowed by the 
experimental data. However, the muon g — 2 [16| (using 
the e + e~ data) restricts mo < 500 GeV at the 2 sigma 
level when tan/3 = 10. 

In conclusion, we have studied the correlation between 
r — > /i7 and B s -B s mixing using SU(5) and SO(10) mod- 
els. The SO(10) GUT models can have larger effects on 
B s -B s mixing compared to the SU(5) boundary condi- 
tions. We find that when Br(r — > /17) is enhanced around 
the current experimental bound, B s -B s mixing may also 
contain large contribution from SUSY. Consequently, the 
phase of Mi2(B s ) is large for a SO(10) model and can be 
tested by measuring CP asymmetries of the B s decay 
modes. It is interesting to note that the phase of B s -B s 
mixing can be large for smaller tan/3 and large scalar 
mass mo at the GUT scale. 
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